-1- 

METHOD AND APPARATUS FOR MULT I -SENSORY 
SPEECH ENHANCEMENT 

BACKGROUND OF THE INVENTION 
The present invention relates to noise 
5 reduction. In particular, the present invention 
relates to removing noise from speech signals. 

A common problem in speech recognition and 
speech transmission is the corruption of the speech 
signal by additive noise. In particular, corruption 
10 due to the speech of another speaker has proven to be 
difficult to detect and/or correct. 

One technique for removing noise attempts 
to model the noise using a set of noisy training 
signals collected under various conditions. These 
15 training signals are received before a test signal 
that is to be decoded or transmitted and are used for 
training purposes only. Although such systems attempt 
to build models that take noise into consideration, 
they are only effective if the noise conditions of 
20 the training signals match the noise conditions of 
the test signals. Because of the large number of 
possible noises and the seemingly infinite 
combinations of noises, it is very difficult to build 
noise models from training signals that can handle 
25 every test condition. 

Another technique for removing noise is to 
estimate the noise in the test signal and then 
subtract it from the noisy speech signal. Typically, 
such systems estimate the noise from previous frames 
30 of the test signal. As such, if the noise is 
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changing over time, the estimate of the noise for the 
current frame will be inaccurate. 

One system of the prior art for estimating 
the noise in a speech signal uses the harmonics of 
5 human speech- The harmonics of human speech produce 
peaks in the frequency spectrum. By identifying 
nulls between these peaks, these systems identify the 
spectrum of the noise. This spectrum is then 
subtracted from the spectrum of the noisy speech 

10 signal to provide a clean speech signal. 

The harmonics of speech have also been used 
in speech coding to reduce the amount of data that 
must be sent when encoding speech for transmission 
across a digital communication path. Such systems 

15 attempt to separate the speech signal into a harmonic 
component and a random component. Each component is 
then encoded separately for transmission. One system 
in particular used a harmonic+noise model in which a 
sum-of-sinusoids model is fit to the speech signal to 

20 perform the decomposition. 

In speech coding, the decomposition is done 
to find a parameterization of the speech signal that 
accurately represents the input noisy speech signal. 
The decomposition has no noise-reduction capability. 

25 Recently, a system has been developed that 

attempts to remove noise by using a combination of an 
alternative sensor, such as a bone conduction 
microphone, and an air conduction microphone. This 
system is trained using three training channels: a 

30 noisy alternative sensor training signal, a noisy air 
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conduction microphone training signal, and a clean 
air conduction microphone training signal. Each of 
the signals is converted into a feature domain. The 
features for the noisy alternative sensor signal and 
5 the noisy air conduction microphone signal are 
combined into a single vector representing a noisy 
signal. The features for the clean air conduction 
microphone signal form a single clean vector. These 
vectors are then used to train a mapping between the 

10 noisy vectors and the clean vectors. Once trained, 
the mappings are applied to a noisy vector formed 
from a combination of a noisy alternative sensor test 
signal and a noisy air conduction microphone test 
signal. This mapping produces a clean signal vector. 

15 This system is less than optimum when the 

noise conditions of the test signals do not match the 
noise conditions of the training signals because the 
mappings are designed for the noise conditions of the 
training signals. 

20 SUMMARY OF THE INVENTION 

A method and system use an alternative 
sensor signal received from a sensor other than an 
air conduction microphone to estimate a clean speech 
value. The clean speech value is estimated without 

25 using a model trained from noisy training data 
collected from an air conduction microphone. Under 
one embodiment, correction vectors are added to a 
vector formed from the alternative sensor signal in 
order to form a filter, which is applied to the air 

30 conductive microphone signal to produce the clean 
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speech estimate. In other embodiments, the pitch of 
a speech signal is determined from the alternative 
sensor signal and is used to decompose an air 
conduction microphone signal. The decomposed signal 
5 is then used to identify a clean signal estimate. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of one computing 
environment in which the present invention may be 
practiced. 

10 FIG. 2 is a block diagram of an alternative 

computing environment in which the present invention 

may be practiced. 

FIG. 3 is a block diagram of a general 

speech processing system of the present invention. 
15 FIG. 4 is a block diagram of a system for 

training noise reduction parameters under one 

embodiment of the present invention. 

FIG. 5 is a flow diagram for training noise 

reduction parameters using the system of FIG. 4. 
20 FIG. 6 is a block diagram of a system for 

identifying an estimate of a clean speech signal from 

a noisy test speech signal under one embodiment of 

the present invention. 

FIG. 7 is a flow diagram of a method for 
25 identifying an estimate of a clean speech signal 

using the system of FIG. 6. 

FIG. 8 is a block diagram of an alternative 

system for identifying an estimate of a clean speech 

signal . 
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FIG. 9 is a block diagram of a second 
alternative system for identifying an estimate of a 
clean speech signal. 

FIG . 10 is a flow diagram of a method for 
5 identifying an estimate of a clean speech signal 
using the system of FIG. 9. 

FIG. 11 is a block diagram of a bone 
conduction microphone. 

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

10 FIG. 1 illustrates an example of a suitable 

computing system environment 100 on which the 
invention may be implemented. The computing system 
environment 100 is only one example of a suitable 
computing environment and is not intended to suggest 

15 any limitation as to the scope of use or 
functionality of the invention. Neither should the 
computing environment 100 be interpreted as having 
any dependency or requirement relating to any one or 
combination of components illustrated in the 

20 exemplary operating environment 100. 

The invention is operational with numerous 
other general purpose or special purpose computing 
system environments or configurations. Examples of 
well-known computing systems, environments, and/or 

25 configurations that may be suitable for use with the 
invention include, but are not limited to, personal 
computers, server computers, hand-held or laptop 
devices, multiprocessor systems, microprocessor-based 
systems, set top boxes, programmable consumer 

30 electronics, network PCs, minicomputers, mainframe 
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computers, telephony systems, distributed computing 
environments that include any of the above systems or 
devices, and the like. 

The invention may be described in the 
5 general context of computer-executable instructions, 
such as program modules, being executed by a 
computer. Generally, program modules include 
routines, programs, objects, components, data 
structures, etc. that perform particular tasks or 

10 implement particular abstract data types. The 
invention is designed to be practiced in distributed 
computing environments where tasks are performed by 
remote processing devices that are linked through a 
communications network. In a distributed computing 

15 environment, program modules are located in both 
local and remote computer storage media including 
memory storage devices. 

With reference to FIG. 1, an exemplary 
system for implementing the invention includes a 

20 general-purpose computing device in the form of a 
computer 110. Components of computer 110 may 

include, but are not limited to, a processing unit 
120, a system memory 130, and a system bus 121 that 
couples various system components including the 

25 system memory to the processing unit 120. The system 
bus 121 may be any of several types of bus structures 
including a memory bus or memory controller, a 
peripheral bus, and a local bus using any of a 
variety of bus architectures. By way of example, and 

30 not limitation, such architectures include Industry 
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Standard Architecture (ISA) bus, Micro Channel 
Architecture (MCA) bus, Enhanced ISA (EISA) bus, 
Video Electronics Standards Association (VESA) local 
bus, and Peripheral Component Interconnect (PCI) bus 
5 also known as Mezzanine bus. 

Computer 110 typically includes a variety 
of computer readable media. Computer readable media 
can be any available media that can be accessed by 
computer 110 and includes both volatile and 

10 nonvolatile media, removable and non-removable media. 
By way of example, and not limitation, computer 
readable media may comprise computer storage media 
and communication media. Computer storage media 
includes both volatile and nonvolatile, removable and 

15 non-removable media implemented in any method or 
technology for storage of information such as 
computer readable instructions, data structures, 
program modules or other data. Computer storage 
media includes, but is not limited to, RAM, ROM, 

20 EEPROM, flash memory or other memory technology, CD- 
ROM, digital versatile disks (DVD) or other optical 
disk storage, magnetic cassettes, magnetic tape, 
magnetic disk storage or other magnetic storage 
devices, or any other medium which can be used to 

25 store the desired information and which can be 
accessed by computer 110. Communication media 

typically embodies computer readable instructions, 
data structures, program modules or other data in a 
modulated data signal such as a carrier wave or other 

30 transport mechanism and includes any information 
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delivery media. The term "modulated data signal" 
means a signal that has one or more of its 
characteristics set or changed in such a manner as to 
encode information in the signal. By way of example, 
5 and not limitation, communication media includes 
wired media such as a wired network or direct-wired 
connection, and wireless media such as acoustic, RF, 
infrared and other wireless media. Combinations of 
any of the above should also be included within the 

10 scope of computer readable media. 

The system memory 130 includes computer 
storage media in the form of volatile and/or 
nonvolatile memory such as read only memory (ROM) 131 
and random access memory (RAM) 132. A basic 

15 input/output system 133 (BIOS) , containing the basic 
routines that help to transfer information between 
elements within computer 110, such as during start- 
up, is typically stored in ROM 131. RAM ' 132 
typically contains data and/or program modules that 

20 are immediately accessible to and/or presently being 
operated on by processing unit 120. By way of 
example, and not limitation, FIG. 1 illustrates 
operating system 134, application programs 135, other 
program modules 136, and program data 137. 

25 The computer 110 may also include other 

removable/non-removable volatile/nonvolatile computer 
storage media. By way of example only, FIG. 1 
illustrates a hard disk drive 141 that reads from or 
writes to non-removable, nonvolatile magnetic media, 

30 a magnetic disk drive 151 that reads from or writes 
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to a removable, nonvolatile magnetic disk 152, and an 
optical disk drive 155 that reads from or writes to a 
removable, nonvolatile optical disk 156 such as a CD 
ROM or other optical media. Other removable/non- 
5 removable, volatile/nonvolatile computer storage 
media that can be used in the exemplary operating 
environment include, but are not limited to, magnetic 
tape cassettes, flash memory cards, digital versatile 
disks, digital video tape, solid state RAM, solid 

10 state ROM, and the like. The hard disk drive 141 is 
typically connected to the system bus 121 through a 
non-removable memory interface such as interface 140, 
and magnetic disk drive 151 and optical disk drive 
155 are typically connected to the system bus 121 by 

15 a removable memory interface, such as interface 150. 

The drives and their associated computer 
storage media discussed above and illustrated in FIG. 
1, provide storage of computer readable instructions, 
data structures, program modules and other data for 

20 the computer 110. In FIG. 1, for example, hard disk 
drive 141 is illustrated as storing operating system 
144, application programs 145, other program modules 
146, and program data 147. Note that these 

components can either be the same as or different 

25 from operating system 134, application programs 135, 
other program modules 136, and program data 137. 
Operating system 144, application programs 145, other 
program modules 14 6, and program data 147 are given 
different numbers here to illustrate that, at a 

30 minimum, they are different copies. 
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A user may enter commands and information 
into the computer 110 through input devices such as a 
keyboard 162, a microphone 163, and a pointing device 
161, such as a mouse, trackball or touch pad. Other 
5 input devices (not shown) may include a joystick, 
game pad, satellite dish, scanner, or the like. 
These and other input devices are often connected to 
the processing unit 120 through a user input 
interface 160 that is coupled to the system bus, but 

10 may be connected by other interface and bus 
structures, such as a parallel port, game port or a 
universal serial bus (USB) . A monitor 191 or other 
type of display device is also connected to the 
system bus 121 via an interface, such as a video 

15 interface 190. In addition to the monitor, computers 
may also include other peripheral output devices such 
as speakers 197 and printer 196, which may be 
connected through an output peripheral interface 195. 

The computer 110 is operated in a networked 

20 environment using logical connections to one or more 
remote computers, such as a remote computer 180. The 
remote computer 180 may be a personal computer, a 
hand-held device, a server, a router, a network PC, a 
peer device or other common network node, and 

25 typically includes many or all of the elements 
described above relative to the computer 110. The 
logical connections depicted in FIG. 1 include a 
local area network (LAN) 171 and a wide area network 
(WAN) 173, but may also include other networks. Such 

30 networking environments are commonplace in offices, 
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enterprise-wide computer networks, intranets and the 
Internet . 

When used in a LAN networking environment, 
the computer 110 is connected to the LAN 171 through 
5 a network interface or adapter 170. When used in a 
WAN networking environment, the computer 110 
typically includes a modem 172 or other means for 
establishing communications over the WAN 173, such as 
the Internet. The modem 172, which may be internal 

10 or external, may be connected to the system bus 121 
via the user input interface 160, or other 
appropriate mechanism. In a networked environment, 
program modules depicted relative to the computer 
110, or portions thereof, may be stored in the remote 

15 memory storage device. By way of example, and not 
limitation, FIG. 1 illustrates remote application 
programs 185 as residing on remote computer 180. It 
will be appreciated that the network connections 
shown are exemplary and other means of establishing a 

20 communications link between the computers may be 
used. 

FIG. 2 is a block diagram of a mobile 
device 200, which is an exemplary computing 
environment. Mobile device 200 includes a 

25 microprocessor 202, memory 204, input/output (I/O) 
components 206, and a communication interface 208 for 
communicating with remote computers or other mobile 
devices. In one embodiment, the afore-mentioned 
components are coupled for communication with one 

30 another over a suitable bus 210. 
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Memory 204 is implemented as non-volatile 
electronic memory such as random access memory (RAM) 
with a battery back-up module (not shown) such that 
information stored in memory 204 is not lost when the 
5 general power to mobile device 200 is shut down. A 
portion of memory 204 is preferably allocated as 
addressable memory for program execution, while 
another portion of memory 204 is preferably used for 
storage, such as to simulate storage on a disk drive. 

10 Memory 204 includes an operating system 

212, application programs 214 as well as an object 
store 216. During operation, operating system 212 is 
preferably executed by processor 202 from memory 204. 
Operating system 212, in one preferred embodiment ,• is 

15 a WINDOWS® CE brand operating system commercially 
available from Microsoft Corporation. Operating 
system 212 is preferably designed for mobile devices, 
and implements database features that can be utilized 
by applications 214 through a set of exposed 

20 application programming interfaces and methods. The 
objects in object store 216 are maintained by 
applications 214 and operating system 212, at least 
partially in response to calls to the exposed 
application programming interfaces and methods. 

25 Communication interface 208 represents 

numerous devices and technologies that allow mobile 
device 200 to send and receive information. The 
devices include wired and wireless modems, satellite 
receivers and broadcast tuners to name a few. Mobile 

30 device 200 can also be directly connected to a 
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computer to exchange data therewith. In such cases, 
communication interface 208 can be an infrared 
transceiver or a serial or parallel communication 
connection, all of which are capable of transmitting 
5 streaming information. 

Input /output components 206 include a 
variety of input devices such as a touch-sensitive 
screen, buttons, rollers, and a microphone as well as 
a variety of output devices including an audio 

10 generator, a vibrating device, and a display. The 
devices listed above are by way of example and need 
not all be present on mobile device 200. In 
addition, other input/output devices may be attached 
to or found with mobile device 200 within the scope 

15 of the present invention. 

FIG. 3 provides a basic block diagram of 
embodiments of the present invention. In FIG. 3, a 
speaker 300 generates a speech signal 302 that is 
detected by an air conduction microphone 30.4 and an 

20 alternative sensor 306. Examples of alternative 
sensors include a throat microphone that measures the 
user's throat vibrations, a bone conduction sensor 
that is located on or adjacent to a facial or skull 
bone of the user (such as the jaw bone) or in the ear 

25 of the user and that senses vibrations of the skull 
and jaw that correspond to speech generated by the 
user. Air conduction microphone 304 is the type of 
microphone that is used commonly to convert audio 
air-waves into electrical signals. 
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Air conduction microphone 304 also receives 
noise 308 generated by one or more noise sources 310. 
Depending on the type of alternative sensor and the 
level of the noise, noise 308 may also be detected by 
5 alternative sensor 306. However, under embodiments 
of the present invention, alternative sensor 306 is 
typically less sensitive to ambient noise than air 
conduction microphone 304. Thus, the alternative 
sensor signal 312 generated by alternative sensor 306 

10 generally includes less noise than air conduction 
microphone signal 314 generated by air conduction 
microphone 304. 

Alternative sensor signal 312 and air 
conduction microphone signal 314 are provided to a 

15 clean signal estimator 316, which estimates a clean 
signal 318. Clean signal estimate 318 is provided to 
a speech process 320. Clean signal estimate 318 may 
either be a filtered time-domain signal or a feature 
domain vector. If clean signal estimate 318 is a 

20 time-domain signal, speech process 320 may take the 
form of a listener, a speech coding system, or a 
speech recognition system. If clean signal estimate 
318 is a feature domain vector, speech process 320 
will typically be a speech recognition system. 

25 The present invention provides several 

methods and systems for estimating clean speech using 
air conduction microphone signal 314 and alternative 
sensor signal 312. One system uses stereo training 
data to train correction vectors for the alternative 

30 sensor signal. When these correction vectors are 
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later added to a test alternative sensor vector, they 
provide an estimate of a clean signal vector. One 
further extension of this system is to first track 
time-varying distortion and then to incorporate this 
5 information into the computation of the correction 
vectors and into the estimation of clean speech. 

A second system provides an interpolation 
between the clean signal estimate generated by the 
correction vectors and an estimate formed by 

10 subtracting an estimate of the current noise in the 
air conduction test signal from the air conduction 
signal. A third system uses the alternative sensor 
signal to estimate the pitch of the speech signal and 
then uses the estimated pitch to identify an estimate 

15 for the clean signal. Each of these systems is 
discussed separately below. 

TRAINING STEREO CORRECTION VECTORS 
FIGS. 4 and 5 provide a block diagram and 
flow diagram for training stereo correction vectors 

20 for the two embodiments of the present invention that 
rely on correction vectors to generate an estimate of 
clean speech. 

The method of identifying correction 
vectors begins in step 500 of FIG. 5, where a "clean" 

25 air conduction microphone signal is converted into a 
sequence of feature vectors. To do this, a speaker 
400 of FIG. 4, speaks into an air conduction 
microphone 410, which converts the audio waves into 
electrical signals. The electrical signals are then 

30 sampled by an analog-to-digital converter 414 to 
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generate a sequence of digital values, which are 
grouped into frames of values by a frame constructor 
416. In one embodiment, A-to-D converter 414 samples 
the analog signal at 16 kHz and 16 bits per sample, 
5 thereby creating 32 kilobytes of speech data per 
second and frame constructor 416 creates a new frame 
every 10 milliseconds that includes 25 milliseconds 
worth of data. 

Each frame of data provided by frame 

10 constructor 416 is converted into a feature vector by 
a feature extractor 418. Under one embodiment, 
feature extractor 418 forms cepstral features. 
Examples of such features include LPC derived 
cepstrum, and Mel-Frequency Cepstrum Coefficients. 

15 Examples of other possible feature extraction modules 
that may be used with the present invention include 
modules for performing Linear Predictive Coding 
(LPC), Perceptive Linear Prediction (PLP) , and 
Auditory model feature extraction. Note that the 

20 invention is not limited to these feature extraction 
modules and that other modules may be used within the 
context of the present invention. 

In step 502 of FIG. 5, an alternative 
sensor signal is converted into feature vectors. 

25 Although the conversion of step 502 is shown as 
occurring after the conversion of step 500, any part 
of the conversion may be performed before, during or 
after step 500 under the present invention. The 
conversion of step 502 is performed through a process 

30 similar to that described above for step 500. 
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In the embodiment of FIG. 4, this process 
begins when alternative sensor 402 detects a physical 
event associated with the production of speech by 
speaker 400 such as bone vibration or facial 
5 movement. As shown in FIG. 11, in one embodiment of a 
bone conduction sensor 1100, a soft elastomer bridge 
1102 is adhered to the diaphragm 1104 of a normal air 
conduction microphone 1106. This soft bridge 1102 
conducts vibrations from skin contact 1108 of the 

10 user directly to the diaphragm 1104 of microphone 
1106. The movement of diaphragm 1104 is converted 
into an electrical signal by a transducer 1110 in 
microphone 1106. Alternative sensor 402 converts the 
physical event into analog electrical signal, which 

15 is sampled by an analog-to-digital converter 404. 
The sampling characteristics for A/D converter 404 
are the same as those described above for A/D 
converter 414. The samples provided by A/D converter 
404 are collected into frames by a frame constructor 

20 406, which acts in a manner similar to frame 
constructor 416. These frames of samples are then 
converted into feature vectors by a feature extractor 
408, which uses the same feature extraction method as 
feature extractor 418. 

25 The feature vectors for the alternative 

sensor signal and the air conductive signal are 
provided to a noise reduction trainer 420 in FIG. 4. 
At step 504 of FIG. 5, noise reduction trainer 420 
groups the feature vectors for the alternative sensor 

30 signal into mixture components. This grouping can be 
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done by grouping similar feature vectors together 
using a maximum likelihood training technique or by 
grouping feature vectors that represent a temporal 
section of the speech signal together. Those skilled 
5 in the art will recognize that other techniques for 
grouping the feature vectors may be used and that the 
two techniques listed above are only provided as 
examples . 

Noise reduction trainer 420 then determines 
10 a correction vector, r s , for each mixture component, 
s, at step 508 of FIG. 5. Under one embodiment, the 
correction vector for each mixture component is 
determined using maximum likelihood criterion. Under 
this technique, the correction vector is calculated 
15 as : 

5>(s 1 6, )(*,-*>,) 
r s = x-i EQ . 1 

Where x t is the value of the air conduction 
vector for frame t and b t is the value of the 

20 alternative sensor vector for frame t. In Equation 
1: 

where p(s) is simply one over the number of mixture 
components and p(b t \s) is modeled as a Gaussian 
25 distribution: 

p{b t | s) = N(b,;fi b9 r b ) EQ.3 
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with the mean and variance T b trained using an 

Expectation Maximization (EM) algorithm where each 
iteration consists of the following steps: 

Y s (t) = p{s\b,) EQ.4 

5 

2>.(0 

EQ.4 is the E-step in the EM algorithm, which uses 
the previously estimated parameters. EQ.5 and EQ.6 

10 are the M-step, which updates the parameters using 
the E-step results. 

The E- and M-steps of the algorithm iterate 
until stable values for the model parameters are 
determined. These parameters are then used to 

15 evaluate equation 1 to form the correction vectors. 
The correction vectors and the model parameters are 
then stored in a noise reduction parameter storage 
422. 

After a correction vector has been 
20 determined for each mixture component at step 508, 
the process of training the noise reduction system of 
the present invention is complete. Once a correction 
vector has been determined for each mixture, the 
vectors may be used in a noise reduction technique of 
25 the present invention. Two separate noise reduction 
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techniques that use the correction vectors are 
discussed below. 

NOISE REDUCTION USING CORRECTION VECTOR 
AND NOISE ESTIMATE 

5 A system and method that reduces noise in a 

noisy speech signal based on correction vectors and a 
noise estimate is shown in the block diagram of FIG. 
6 and the flow diagram of FIG. 7 , respectively. 

At step 700, an audio test signal detected 

10 by an air conduction microphone 604 is converted into 
feature vectors. The audio test signal received by 
microphone 604 includes speech from a speaker 600 and 
additive noise from one or more noise sources 602. 
The audio test signal detected by microphone 604 is 

15 converted into an electrical signal that is provided 
to analog-to-digital converter 606. 

A-to-D converter 606 converts the analog 
signal from microphone 604 into a series of digital 
values. In several embodiments, A-to-D converter 606 

20 samples the analog signal at 16 kHz and 16 bits per 
sample, thereby creating 32 kilobytes of speech data 
per second. These digital values are provided to a 
frame constructor 607, which, in one embodiment, 
groups the values into 25 millisecond frames that 

25 start 10 milliseconds apart. 

The frames of data created by frame 
constructor 607 are provided to feature extractor 
610, which extracts a feature from each frame. Under 
one embodiment, this feature extractor is different 

30 from feature extractors 408 and 418 that were used to 
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train the correction vectors. In particular, under 
this embodiment, feature extractor 610 produces power 
spectrum values instead of cepstral values. The 
extracted features are provided to a clean signal 
5 estimator 622, a speech detection unit 626 and a 
noise model trainer 624. 

At step 702, a physical event, such as bone 
vibration or facial movement, associated with the 
production of speech by speaker 600 is converted into 

10 a feature vector. Although shown as a separate step 
in FIG. 7, those skilled in the art will recognize 
that portions of this step may be done at the same 
time as step 700. During step 702, the physical 
event is detected by alternative sensor 614. 

15 Alternative sensor 614 generates an analog electrical 
signal based on the physical events. This analog 
signal is converted into a digital signal by analog- 
to-digital converter 616 and the resulting digital 
samples are grouped into frames by frame constructor 

20 617. Under one embodiment, analog-to-digital 

converter 616 and frame constructor 617 operate in a 
manner similar to analog-to-digital converter 606 and 
frame constructor 607. 

The frames of digital values are provided 

25 to a feature extractor 620, which uses the same 
feature extraction technique that was used to train 
the correction vectors. As mentioned above, examples 
of such feature extraction modules include modules 
for performing Linear Predictive Coding (LPC) , LPC 

30 derived cepstrum, Perceptive Linear Prediction (PLP), 
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Auditory model feature extraction, and Mel-Frequency 
Cepstrum Coefficients (MFCC) feature extraction. In 
many embodiments, however, feature extraction 
techniques that produce cepstral features are used. 
5 The feature extraction module produces a 

stream of feature vectors that are each associated 
with a separate frame of the speech signal. This 
stream of feature vectors is provided to clean signal 
estimator 622. 

10 The frames of values from frame constructor 

617 are also provided to a feature extractor 621, 
which in one embodiment extracts the energy of each 
frame. The energy value for each frame is provided 
to a speech detection unit 626. 

15 At step 704, speech detection unit 626 uses 

the energy feature of the alternative sensor signal 
to determine when speech is likely present. This 
information is passed to noise model trainer 624, 
which attempts to model the noise during periods when 

20 there is no speech at step 706. 

Under one embodiment, speech detection unit 
626 first searches the sequence of frame energy 
values to find a peak in the energy. It then 
searches for a valley after the peak. The energy of 

25 this valley is referred to as an energy separator, d. 
To determine if a frame contains speech, the ratio, 
k, of the energy of the frame, e, over the energy 
separator, d, is then determined as: k=e/d. A speech 
confidence, q, for the frame is then determined as: 
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where a defines the transition between two states and 
in one implementation is set to 2. Finally, we use 
the average confidence value of its 5 neighboring 
5 frames (including itself) as the final confidence 
value for this frame. 



value is used to determine if speech is present such 
that if the confidence value exceeds the threshold, 
10 the frame is considered to contain speech and if the 
confidence value does not exceed the threshold, the 
frame is considered to contain non-speech. Under one 
embodiment, a threshold value of 0.1 is used. 



15 speech detection unit 626, noise model trainer 624 
updates a noise model 625 at step 706. Under one 
embodiment, noise model 625 is a Gaussian model that 
has a mean fd n and a variance I„ . This model is based 
on a moving window of the most recent frames of non- 
20 speech. Techniques for determining the mean and 
variance from the non-speech frames in the window are 
well known in the art. 



parameter storage 422 and noise model 625 are 
25 provided to clean signal estimator 622 with the 
feature vectors, b, for the alternative sensor and 



Under one embodiment, a fixed threshold 



For each non-speech frame detected by 



Correction vectors and model parameters in 
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the feature vectors, S y , for the noisy air conduction 

microphone signal. At step 708, clean signal 

estimator 622 estimates an initial value for the 
clean speech signal based on the alternative sensor 
5 feature vector, the correction vectors, and the model 
parameters for the alternative sensor. In 
particular, the alternative sensor estimate of the 
clean signal is calculated as: 

x = b + ^p(s\b)r s EQ.8 

s 

10 where x is the clean signal estimate in the cepstral 
domain, b is the alternative sensor feature vector, 
p{s\b) is determined using equation 2 above, and r s is 

the correction vector for mixture component s. Thus, 
the estimate of the clean signal in Equation 8 is 
15 formed by adding the alternative sensor feature 
vector to a weighted sum of correction vectors where 
the weights are based on the probability of a mixture 
component given the alternative sensor feature 
vector. 

20 At step 710, the initial alternative sensor 

clean speech ' estimate is refined by combining it with 
a clean speech estimate that is formed from the noisy 
air conduction microphone vector and the noise model. 
This results in a refined clean speech estimate 628. 

25 In order to combine the cepstral value of the initial 
clean signal estimate with the power spectrum feature 
vector of the noisy air conduction microphone, the 
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cepstral value is converted to the power spectrum 
domain using: 

S xlb =e c ' ] " EQ.9 

where C" 1 is an inverse discrete cosine transform and 

5 S xlb is the power spectrum estimate of the clean signal 

based on the alternative sensor. 

Once the initial clean signal estimate from 
the alternative sensor has been placed in the power 
spectrum domain, it can be combined with the noisy 
10 air conduction microphone vector and the noise model 
as : 

s, = (s; 1 + z J r' rc' (s y -m„)+ r;;^ ] eq . i o 

where S x is the refined clean signal estimate in the 
power spectrum domain, S y is the noisy air conduction 
15 microphone feature vector, (// M> S M ) are the mean and 
covariance of the prior noise model (see 624), S x[b is 
the initial clean signal estimate based on the 
alternative sensor, and is the covariance matrix 

of the conditional probability distribution for the 
20 clean speech given the alternative sensor's 
measurement. can be computed as follows. 

LetJdenote the Jacobian of the function on the right 
hand side of equation 9. Let S be the covariance 

matrix of x . Then the covariance of S xlb is 
25 l Mb =JU T EQ. 11 
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In a simplified embodiment, we rewrite 
EQ.10 as the following equation: 

S x =a(f)(S y -M„) + (l-a(f))S xlb EQ. 12 

5 where cc{f) is a function of both the time and the 
frequency band. Since the alternative sensor that we 
are currently using has the bandwidth up to 3KHz, we 
choose &(f) to be 0 for the frequency band below 
3KHz. Basically, we trust the initial clean signal 

10 estimate from the alternative sensor for low 
frequency bands. For high frequency bands, the 
initial clean signal estimate from the alterative 
sensor is not so reliable. Intuitively, when the 
noise is small for a frequency band at the current 

15 frame, we would like to choose a large oc(f) so that 
we use more information from the air conduction 
microphone for this frequency band. Otherwise, we 
would like to use more information from the 
alternative sensor by choosing a small oc{f) . In one 

20 embodiment, we use the energy of the initial clean 
signal estimate from the alternative sensor to 
determine the noise level for each frequency band. 
Let E(f) denote the energy for frequency band / . Let 

M-Max f E{f). oc(f) , as a function of /, is defined as 

25 follows: 



«(/)= 
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f-3K 



: f>4K 

M 



a(4K) : 3K<f<4K EQ. 13 



1/: 

0 : /<3£ 



where we use a linear interpolation to transition 
from 3K to 4K to ensure the smoothness of a(f) . 

The refined clean signal estimate in the 
5 power spectrum domain may be used to construct a 
Wiener filter to filter the noisy air conduction 
microphone signal. In particular, the Wiener filter, 
H, is set such that: 

H = — EQ.14 
Sy 

10 This filter can then be applied against the 

time domain noisy air conduction microphone signal to 
produce a noise-reduced or clean time-domain signal. 
The noise-reduced signal can be provided to a 
listener or applied to a speech recognizer. 

15 Note that Equation 12 provides a refined 

clean signal estimate that is the weighted sum of two 
factors, one of which is a clean signal estimate from 
an alternative sensor. This weighted sum can be 
extended to include additional factors for additional 

20 alternative sensors. Thus, more than one alternate 
sensor may be used to generate independent estimates 
of the clean signal. These multiple estimates can 
then be combined using equation 12. 
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NOISE REDUCTION USING CORRECTION VECTOR 
WITHOUT NOISE ESTIMATE 

FIG. 8 provides a block diagram of an 
5 alternative system for estimating a clean speech 
value under the present invention. The system of 
FIG. 8 is similar to the system of FIG. 6 except that 
the estimate of the clean speech value is formed 
without the need for an air conduction microphone or 

10 a noise model. 

In FIG. 8, a physical event associated with 
a speaker 800 producing speech is converted into a 
feature vector by alternative sensor 802, analog-to- 
digital converter 804, frame constructor 806 and 

15 feature extractor 808, in a manner similar to that 
discussed above for alternative sensor 614, analog- 
to-digital converter 616, frame constructor 617 and 
feature extractor 618 of FIG. 6. The feature vectors 
from feature extractor 808 and the noise reduction 

20 parameters 422 are provided to a clean signal 
estimator 810, which determines an estimate of a 

clean signal value 812, S xlb , using equations 8 and 9 

above . 

The clean signal estimate, S x{b , in the power 

25 spectrum domain may be used* to construct a Wiener 
filter to filter a noisy air conduction microphone 
signal. In particular, the Wiener filter, H, is set 
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such that: 

H = ^ EQ.15 

s, 

This filter can then be applied against the 
time domain noisy air conduction microphone signal to 
5 produce a noise-reduced or clean signal. The noise- 
reduced signal can be provided to a listener or 
applied to a speech recognizer . 

Alternatively, the clean signal estimate in 
the cepstral domain, x , which is calculated in 
10 Equation 8, may be applied directly to a speech 
recognition system. 

NOISE REDUCTION USING PITCH TRACKING 

An alternative technique for generating 
15 estimates of a clean speech signal is shown in the 
block diagram of FIG. 9 and the flow diagram of FIG. 
10. In particular, the embodiment of FIGS. 9 and 10 
determine a clean speech estimate by identifying a 
pitch for the speech signal using an alternative 
20 sensor and then using the pitch to decompose a noisy 
air conduction microphone signal into a harmonic 
component and a random component. Thus, the noisy 

signal is represented as: 

y = y*+y, eq. 16 

25 where y is the noisy signal, y h is the harmonic 
component, and y r is the random component. A weighted 
sum of the harmonic component and the random 
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component are used to form a noise-reduced feature 
vector representing a noise-reduced speech signal. 

Under one embodiment , the harmonic 
component is modeled as a sum of harmonically-related 
5 sinusoids such that: 

K 

y h = y £ j a k cos(ka) Q t) + b k s\n{kco 0 t) EQ . 17 

k=\ 

where co 0 is the fundamental or pitch frequency and K 
is the total number of harmonics in the signal. 

Thus, to identify the harmonic component, 
10 an estimate of the pitch frequency and the amplitude 
parameters {a x a 2 ...a k b x b 2 ...b k } must be determined. 

At step 1000, a noisy speech signal is 
collected and converted into digital samples. To do 
this, an air conduction microphone 904 converts audio 

15 waves from a speaker 900 and one or more additive 
noise sources 902 into electrical signals. The 
electrical signals are then sampled by an analog-to- 
digital converter 906 to generate a sequence of 
digital values. In one embodiment, A-to-D converter 

20 906 samples the analog signal at 16 kHz and 16 bits 
per sample, thereby creating 32 kilobytes of speech 
data per second. At step 1002, the digital samples 
are grouped into frames by a frame constructor 908. 
Under one embodiment, frame constructor 908 creates a 

25 new frame every 10 milliseconds that includes 25 
milliseconds worth of data. 

At step 1004, a physical event associated 
with the production of speech is detected by 
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alternative sensor 944. In this embodiment, an 
alternative sensor that is able to detect harmonic 
components, such as a bone conduction sensor, is best 
suited to be used as alternative sensor 944. Note 
5 that although step 1004 is shown as being separate 
from step 1000, those skilled in the art will 
recognize that these steps may be performed at the 
same time. The analog signal generated by 

alternative sensor 944 is converted into digital 

10 samples by an analog-to-digital converter 94 6. The 
digital samples are then grouped into frames by a 
frame constructer 948 at step 1006. 

At step 1008, the frames of the alternative 
sensor signal are used by a pitch tracker 950 to 

15 identify the pitch or fundamental frequency of the 
speech. 

An estimate for the pitch frequency can be 
determined using any number of available pitch 
tracking systems. Under many of these systems, 

20 candidate pitches are used to identify possible 
spacing between the centers of segments of the 
alternative sensor signal. For each candidate pitch, 
a correlation is determined between successive 
segments of speech. In general, the candidate pitch 

25 that provides the best correlation will be the pitch 
frequency of the frame. In some systems, additional 
information is used to refine the pitch selection 
such as the energy of the signal and/or an expected 
pitch track. 
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Given an estimate of the pitch from pitch 
tracker 950, the air conduction signal vector can be 
decomposed into a harmonic component and a random 
component at step 1010. To do so, equation 17 is 
rewritten as: 

y = Ab EQ. 18 

where y is a vector of N samples of the noisy speech 
signal, A is an Nx2K matrix given by: 

A = [A C0S A sin ] EQ. 19 

with elements 

A cos (*,0 = cos(*6> 0 0 A sin (kj) = sm(ko) 0 t) EQ. 20 

and b is a 2Kx\ vector given by: 



b 7 ' =[a l a 2 ...a k b l b 2 ...b k ] 



EQ. 21 



Then, the least-squares solution for the amplitude 
coefficients is : 

b = (A r Ar , A r y 



Using b, 
component of the 
determined as: 



EQ. 22 

an estimate for the harmonic 
noisy speech signal can be 



y h =Ab EQ. 23 

An estimate of the random component is then 
calculated as: 

y,=y-y* eq. 24 

Thus, using equations 18-24 above, harmonic 
decompose unit 910 is able to produce a vector of 
harmonic component samples 912, y h , and a vector of 
random component samples 914, y r . 
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After the samples of the frame have been 
decomposed into harmonic and random samples, a 
scaling parameter or weight is determined for the 
harmonic component at step 1012. This scaling 
5 parameter is used as part of a calculation of a 
noise-reduced speech signal as discussed further 
below. Under one embodiment, the scaling parameter 
is calculated as: 

ah =ZM>L EQ. 25 

10 where a h is the scaling parameter, y h (i) is the ith 
sample in the vector of harmonic component samples 
and y(i) is the ith sample of the noisy speech signal 
for this frame. In Equation 25, the numerator is the 
sum of the energy of each sample of the harmonic 

15 component and the denominator is the sum of the 
energy of each sample of the noisy speech signal. 
Thus, the scaling parameter is the ratio of the 
harmonic energy of the frame to the total energy of 
the frame. 

20 In alternative embodiments, the scaling 

parameter is set using a probabilistic voiced- 
unvoiced detection unit. Such units provide the 
probability that a particular frame of speech is 
voiced, meaning that the vocal cords resonate during 

25 the frame, rather than unvoiced. The probability 
that the frame is from a voiced region of speech can 
be used directly as the scaling parameter. 
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After the scaling parameter has been 
determined or while it is being determined, the Mel 
spectra for the vector of harmonic component samples 
and the vector of random component samples are 
5 determined at step 1014. This involves passing each 
vector of samples through a Discrete Fourier 
Transform (DFT) 918 to produce a vector of harmonic 
component frequency values 922 and a vector of random 
component frequency values 920. The power spectra 

10 represented by the vectors of frequency values are 
then smoothed by a Mel weighting unit 924 using a 
series of triangular weighting functions applied 
along, the Mel scale. This results in a harmonic 
component Mel spectral vector 928, , and a random 

15 component Mel spectral vector 926, Y r . 

At step 1016, the Mel spectra for the 
harmonic component and the random component are 
combined as a weighted sum to form an estimate of a 
noise-reduced Mel spectrum. This step is performed 

20 by weighted sum calculator 930 using the scaling 
factor determined above in the following equation: 

X(t) = a h (t)Y h (t) + a r Y r (0 EQ. 26 

A 

where X(7) is the estimate of the noise-reduced Mel 
spectrum, Y h (t) is the harmonic component Mel 
25 spectrum, Y r (t) is the random component Mel spectrum, 
a h (t) is the scaling factor determined above, a r is a 
fixed scaling factor for the random component that in 
one embodiment is set equal to .1, and the time index 
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t is used to emphasize that the scaling factor for 
the harmonic component is determined for each frame 
while the scaling factor for the random component 
remains fixed. Note that in other embodiments, the 
5 scaling factor for the random component may be 
determined for each frame. 

After the noise-reduced Mel spectrum has 
been calculated at step 1016, the log 932 of the Mel 
spectrum is determined and then is applied to a 

10 Discrete Cosine Transform 934 at step 1018. This 
produces a Mel Frequency Cepstral Coefficient (MFCC) 
feature vector 936 that represents a noise-reduced 
speech signal. 

A separate noise-reduced MFCC feature 

15 vector is produced for each frame of the noisy 
signal. These feature vectors may be used for any 
desired purpose including speech enhancement and 
speech recognition. For speech enhancement, the MFCC 
feature vectors can be converted into the power 

20 spectrum domain and can be used with the noisy air 
conduction signal to form a Weiner filter. 

Although the present invention has been 
described with reference to particular embodiments, 
workers skilled in the art will recognize that 

25 changes may be made in form and detail without 
departing from the spirit and scope of the invention. 



